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Abstract 

We investigate the specific heat at constant volume Cy in the BCS (Bardeen-Cooper-Schrieffer)-BEC 
(Bose-Einstein condensation) crossover regime of an ultracold Eermi gas above the supertluid phase tran¬ 
sition temperature T^. Within the framework of the strong-coupling theory developed by Nozieres and 
Schmitt-Rink, we show that this thermodynamic quantity is sensitive to the stability of preformed Cooper 
pairs. That is, while Cy(r > Tc) in the unitary regime is remarkably enhanced by metastable preformed 
Cooper pairs or pairing fluctuations, it is well described by that of an ideal Bose gas of long-lived stable 
molecules in the strong-coupling BEC regime. Using these results, we identify the region where the system 
may be viewed as an almost ideal Bose gas of stable pairs, as well as the pseudogap regime where the system 
is dominated by metastable preformed Cooper pairs, in the phase diagram of an ultracold Eermi gas with 
respect to the strength of a pairing interaction and the temperature. We also show that the calculated specific 
heat agrees with the recent experiment on a ^Ei unitary Eermi gas. Since the formation of preformed Cooper 
pairs is a crucial key in the BCS-BEC crossover phenomenon, our results would be helpful in considering 
how fluctuating preformed Cooper pairs appear in a Eermi gas, to eventually become stable, as one passes 
through the BCS-BEC crossover region. 

PACS numbers: 03.75.Ss, 03.75.-b, 03.70.-Hk 
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I. INTRODUCTION 


The formation of preformed Cooper pairs is a erueial key in eonsidering the BCS (Bardeen- 
Cooper-Sehrieffer)-BEC (Bose-Einstein eondensation) erossover phenomenon [|ll-ll ill, where the 
eharaeter of a Eermi superfluid eontinuously ehanges from the weak-eoupling BCS type to the 


BEC of tightly bound moleeules with inereasing the strength o 


the realization of this erossover phenomenon in "^°K [|12l1 and ^Ei [|13l - ll5l1 Eermi gases, by 


a pairing interaetion. Sinee 


using 


a Eeshbaeh resonaneelllbl-llSll. it has extensively been diseussed both theoretieallv 1 9l - |29ll and 


experimentally 1 I 30 I 442 I 1 how metastable preformed Cooper pairs (that are also referred to in the 
literature as pairing fluetuations) appear in a Eermi gas, to eventually beeome long-lived stable 
pairs, as one passes through the BCS-BEC erossover region above the superfluid phase transition 
temperature T^. Sinee the BCS-BEC erossover is also eonsidered as a erueial key in the fields of 
high-Tc euprates llOLlij], as well as iron based supereonduetors ll44ll . elueidating strong-eoupling 
properties of an ultraeold Eermi gas in this regime would also eontribute to the understanding of 
these strongly eorrleated eleetron systems. 

Although there is no elear phase boundary between the weak-eoupling BCS regime and the 
strong-eoupling BEC regime, it is still an interesting problem to physieally identify the region 
where preformed Cooper pairs dominate over system properties. In this regard, we note that, 
when preformed Cooper pairs appear in a normal Eermi gas, single-partiele Eermi exeitations are 
expeeted to have a gap-like strueture, refleeting their finite dissoeiation energy. This so-ealled 


preformed pair seenario has been diseussed in high-Tc euprates [l45ll . as a possible meehanism of 


the pseudogapped density of states observed in the under-doped regime of this strongly eorrelated 
eleetron svstem [|46l - l49n . In this field, the temperature T* below whieh the pseudogap appears in 
the density of states p{oj) is ealled the pseudogap temperature. Although T* is not aeeompanied 
by any phase transition, it is a useful eharaeteristie temperature to distinguish between the normal 
Eermi liquid regime and the pseudogap regime in the phase diagram of high-Tc euprates. 


In high-Tc euprates, the validity of the preformed pair seenario is still in debate[l45 


5QM52D, 


due to the eomplexity of this system. In eontrast, sinee an ultraeold Eermi gas in the BCS-BEC 
erossover region is simply dominated by strong pairing fluetuations, the preformed pair seenario 
is validated. Indeed, it has been pointed out bbL 13611 that the baek-bending eurve of the single- 
partiele dispersion observed by a reeent photoemission-type experiment on a unitaryFermi 
gas Odll may be a signature of the pseudogap phenomenon. It has also been shown OTL 12911 that 
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the anomalous suppression of the uniform spin susoeptibility;t's observed in a ^Li Fermi gas above 
FcldOfl can be explained as an effect of fluctuating spin-singlet preformed pairs. At present, al¬ 
though the pseudogap temperature T* has not experimentally been determined in an ultracold 
Fermi gas, the existence of this characteristic temperature has theoretically been predicted from 
the calculated density of states p(mf [l24n . As in the case of high-r^ cuprates, T* in an ultracold 
Fermi gas is not accompanied by any phase transition. However, it physically gives the boundary 
between a normal Fermi gas regime and the pseudogap regime, being dominated by fluctuating 
preformed Cooper pairs, in the phase diagram of an ultracold Fermi gas above T^. 

A similar characteristic temperature called the spin-gap temperature, has also been 
predicted [1291]. is determined as the temperature below which the spin susceptibility ;^s in the 
normal state is anomalously suppressed by spin-singlet preformed Cooper pairs. Although is 
not exactly the same as the pseudogap temperature T* they have essentially the same background 
physics, and thus also has the meaning of the boundary between the normal Fermi gas reg ime 


and the preformed-pair regime. We brie% note that this so-called spin-gap phenomenon [SSO has 


also been discussed in high-Fc cuprates [|54. 


550. 


Although T* and conveniently give the boundary around which metastable preformed 
Cooper pairs start to dominate over the system in the BCS-BEC crossover region, they do not 
have any information about where these fluctuating preformed pairs become stable in the strong¬ 
coupling BEC regime. In determining this second boundary, however, the low-energy density of 
states p(m ~ 0) (which gives the pseudogap temperature T*), as well as the spin susceptibility 
(which gives the spin gap temperature T^), are not useful, because both quantities almost vanish 
deep inside BEC regime where most Fermi atoms form spin-singlet bound molecules with a large 
binding energy. 

In this regard, the specific heat at constant volume Cy is promising, because it is finite in the 
whole BCS-BEC crossover region. In addition, Cy is sensitive to the quantum statistics of particles 
in the system, in the sense that, while Cy exhibits a linear-temperature dependence in a Fermi gas, 
it increases with decreasing the temperature in an ideal Bose gas. Furthermore, the specific heat 
has recently become accessible in cold Fermi gas physics 114ill . Thus, the above-mentioned second 
boundary may be determined by using this thermodynamic quantity. 

The purpose of this paper is to theoretically confirm this expectation, to distinguish between 
the pseudogap regime, which is dominated by metastable preformed Cooper pairs, or pairing fluc¬ 
tuations, and the region that can be viewed as a gas of long-lived stable pairs, in the phase diagram 
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of an ultracold Fermi gas. Including pairing fluctuations above within the framework of the 
strong-coupling theory developed by Nozieres and Schmitt-RinkyO, we show that the temperature 
dependence of the specific heat is very different in between the BCS-BEC crossover region and 
the strong-coupling BEC regime. Using this difference, we determine a characteristic temperature 
t which conveniently gives the boundary between the pseudogap regime and the region of stable 
pairs. In addition, the specific heat is also shown to be able to determine the boundary between the 
normal Fermi gas regime and the pseudogap regime. The characteristic temperature f giving the 
latter boundary is found to be consistent with the previous pseudogap temperature T*, as well as 
the spin-gap temperature that are, respectively, obtained from the density of states p(cu) and the 
spin susceptibility's. We also show that that our result on Cy agrees with the recent experiment 
on a ^Li unitary Fermi gas 114 ill . We briefly note that the specific heat in a unitary Fermi gas has 
also been discussed within a T-matrix approximation1. as well as within the combined NSR 
theory with local density approximation yOD. 

This paper is organized as follows. In Sec. II, we explain our strong coupling formalism used 
to calculate the specific heat at constant volume Cy in the BCS-BEC crossover region above T^. In 
Sec. Ill, we show our numerical results on Cy over the entire BCS-BEC crossover region. Here, we 
explain how to determine t and f from the temperature dependence of Cy. Using these character¬ 
istic temperatures, we identify the region where the system is dominated by fluctuating metastable 
preformed Cooper pairs, as well as the region where the system is dominated by long-lived stable 
molecules, in the phase diagram of an ultracold Fermi gas with respect to the interaction strength 
and the temperature. Throughout this paper, we set /z = = U and the system volume V is taken 

to be unity, for simplicity. 


II. FORMULATION 


We consider a two-component uniform Fermi gas in the normal state, described by the BCS 
Hamiltonian, 




P^P.D-^P.O- 


u 


p,<T 


z 

p,p'.q 


^p+(j/2,T^-p+<?/2,J,^- 


■p'+q/2,lCp'+q/2,'\ 5 


( 1 ) 


where Cp^o- is the annihilation operator of a Fermi atom with pseudospin cr ='f, describing two 
atomic hyperfine states, ^p = Sp- n = I (2m) -/i is the kinetic energy, measured from the Fermi 
chemical potential p (where m is an atomic mass). -U is an 5-wave pairing interaction, which we 
treat as a tunable parameter. As usual, we measure the interaction strength in terms of the 5-wave 
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FIG. 1: Feynman diagrams describing the correction term Gnsr to the thermodynamic potential Q. in 
the strong-coupling NSR theoryyD. The solid line is the bare single-particle thermal Green’s function 
G~^{p, i(jJn) - i(^n - (where is the fermion Matsubara frequency), and the dashed line is the attractive 
pairing interaction -U. 


scattering length whieh is related to the pairing interaetion -U as, 


Anas 

m 


U 




1 ’ 
2sp 


( 2 ) 


where is a momentum eutoff. 

We include pairing fluctuations within the ordinary NSR theorylSfl. In this BCS-BEC erossover 
theory, the thermodynamie potential Q = Qq + ^^nsr eonsists of the non-interaeting part. 


Qo = -2r^ln[l+e-^-/^], (3) 

p 

and the fluctuation correction Qnsr> the latter of whieh is diagrammatically given in Fig. [T] The 
summation of these diagrams gives 


f^NSR = -T^ lnr(qr, ivn). 




where v„ is the boson Matsubara frequeney, and 


r{q, iVn) 


m 


Ana, 


-I- 




is the NSR partiele-partiele seattering matrix. Here, 

1 “ fi^p+qll) - f{^-p+qll) 


Uiq, iVn) = - ^ 


^p+q/2 ^-p+ql2 


(4) 


(5) 


( 6 ) 
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is the lowest-order pair-correlation function, describing fluctuations in the Cooper channel (where 
f{x) is the Fermi distribution function). 

We calculate the specific heat at constant volume Cy from the formula, 


Cv 


{-] 

\dTk, 


(7) 


Here, the internal energy E is obtained from Q = Qq + ^nsr via the Legendre transformation. 


E = Q-r 






d d 

T ^n(p, iv„) -h iu—U(q, IVn) 
oT oju 


( 8 ) 


= 2^£p/(^p)-r^r(g,/yj 

P qJVn 

The Fermi chemical potential /i in Eq. ([8]) is determined from the equation for the total number N 
of Fermi atoms, given by 


!dQ. 


2 Yj iVn) 


Q,tVn 


dn 




NSRj 


(9) 


where A^p and A^nsr represent the non-interacting part and the NSR strong-coupling corrections, 
respectively. In this paper, we numerically evaluate Eq. (|7]) from the internal energies E{T) and 


E{T + 5T). 

In the NSR theory [|3l-10], the equation for the sup erfluid phase transition temperature is con¬ 
veniently obtained from the Thouless criterion [jS^, stating that the superfluid instability occurs 
when the particle-particle scattering matrix F(q, zv„) in Eq. ([5]) has a pole at q = v„ = 0. The 
resulting Tc-equation has the same form as the mean-field BCS gap equation at Tc, as 


1 = 


Ana, 


m 


E 


Itanh^ 




IT 2s, 


Pi 


( 10 ) 


Eollowing the standard NSR approachyl-l^, we numerically solve the Tc-equation (fTOl) . together 
with the number equation (HI), to self-consistently determine and ju{Tc) in the BCS-BEC 
crossover region. Above T^, we only deal with the number equation (HI) to determine which 
is used to evaluate the specific heat Cy. 
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FIG. 2: (Color online) (a) Calculated specific heat Cy at Tc, as a function of the interaction strength mea¬ 
sured in terms of the inverse scattering length (k^asy^, normalized by the Fermi momentum kj;. We also 
plot the specific heat C^{T - Ty in a free Fermi gas, as well as the specific heat C® in Eq. ^Y2\ . (b) 
The number of stable molecules N'q at T^. N^c is the contribution from scattering states at T^- (c) Fermi 
chemical potential /i(r = Tc), normalized by the Fermi energy gp- The dashed line shows -E^^^/2, where 
Ebind = 1 1(mol) is the binding energy of a two-body bound molecule. 

III. SPECIFIC HEAT IN THE BCS-BEC CROSSOVER REGION ABOVE Tc 

Figure 121; a) shows the speeifie heat at eonstant volume Cy in the BCS-BEC erossover regime 
of an ultraeold Fermi gas at T = Tc. As expeeted, Cy in the weak-eoupling BCS regime 
((kpai)"' < -1, where kp is the Fermi momentum), as well as that in the strong-coupling BEC 
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FIG. 3: (Color online) Calculated specific heat Cy, as a function of the temperature in the BCS-BEC 
crossover regime of an ultracold Fermi gas above Tc- The dashed line shows the result at T^. 


regime {{k^as) * > 1) are, respeetively, well deseribed by the speeifie heat in a free Fermi gas 115711 . 


n 


C'yiT «T^) = -\ — \N 


( 11 ) 


(where Tp is the Fermi temperature), and the speeifie heat in an ideal Bose gas with At/2 molecules 


at the BEC phase transition temperature Tbec = 0.218rFll3l-l5, 


m 


CUT = Tbec) 


0.963V, 

4 ^ /'(3/2) 


( 12 ) 


where /^(3/2) = 2.612 and ^(5/2) = 1.341 are zeta functions. Although Cy continuously changes 
from Cy to Cy in the BCS-BEC crossover, it experiences anomalous enhancement in the unitary 
regime ((kpas)~^ ~ 0), as seen in Fig. WiU)- 

This remarkable enhancement of Cy originates from the suppression of the entropy 5 = In VF 
by the appearance of preformed Cooper pairs near T^. Since the number of possible micro-states 
VF in a gas of bound molecules with nearly zero center of mass momentum is smaller than VF in 
a simple unbound Fermi gas, the gradual formation of preformed Cooper pairs with decreasing 
temperature nearing suppresses the entropy S. When this suppression is more remarkable at 
lower temperatures, the thermodynamic formula. 


Cv 


-d) 


V,N 


(13) 
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FIG. 4: (Color online) Specific heat Cy in the EEC side, 0.2 < (kpas) ^ < 0.6. 


immediately gives the enhaneement of Cy. In the unitary regime, sueh preformed-pair formation 
oeeurs near Tc, so that the amplifieation of Cy is also restrieted to the region near T^, as shown 
Fig. [31 

To see whether these preformed Cooper pairs are stable or fluetuating, it is eonvenient to divide 
the NSR eorreetion term A^nsr in the number equation ® into the sum of twiee the number Nb of 
stable moleeules and the so-ealled seattering part A^scH] involving eontribution from fluetuating 
metastable preformed pairsj^ 59]. The resulting expression for the total number N of Fermi atoms 


has the form. 


N 


iV" + 2iVB+iVsc. 


(14) 


(For detailed expressions for Nb and see the Appendix.) Then, we And in Fig. [Itb) that 
there are no stable preformed pairs {Nb = 0) in the unitary regime where Cy is remarkably 
amplifled libOn . whieh means that this enhaneement is due to the inerease of metastable preformed 
Cooper pairs or pairing fluetuations. 

Figurel^b) also indieates that the strong-eoupling EEC regime ^'O.S) is dominated by 

long-lived stable moleeules {Nb - N/2). This naturally explains why Cy(r = T^) in this regime 
is well deseribed by C® in Eq. (IT^ . As shown in Fig. l2tc), this result is also eonsistent with the 
well-known result for the Fermi ehemieal potential n that it beeomes negative in the EEC regime. 


and the magnitude |/i| approaehes 


ha 


state in the strong-eoupling limityj-lb 


the binding energy , = 1 /{mal) of a two-body bound 


m. 


As seen in Eig. |3l the amplifieation of Cy in the unitary regime near Tc disappears, as one moves 
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FIG. 5: (Color online) (al)-(cl) Specific heat Cv(r > Tc) in the strong-coupling EEC regime > 1)- 

Cy is the specific heat in an ideal gas of N12 bosons with a molecular mass M - 2m. (a2)-(c2) The number 
A^b of stable pairs. The characteristic temperature f is given as the temperature at which Cy takes a minimal 
value in the EEC regime. 


to the EEC side {{kyasT^ ^'O). To see this more elearly, we summarize in Fig. |4]the temperature 
dependenee of Cy slightly in the EEC side. Noting that long-lived stable moleeules appear when 
(kpasT^ >0.3 (see Fig. Sb)), we expeet that the thermal dissoeiation of these moleeules (with a 
relatively small binding energy Ebind) is responsible for the temperature dependenee of Cy near Tc 
in this regime. Indeed, simply taking into aeeount this effeet by dealing with a two-level system 
with energy a> = 0 and to = Ebind, one has 


Cy = 



^bind 

^ 7 ^ 


(15) 


whieh monotonieally inereases with inereasing the temperature when T < £'bind/2. This is just the 
behavior of Cy(r >Tc) shown in Fig. |4]when (kpa^)"^ 2* 0.3, indicating that the increase of Cy 
with increasing the temperature near in this region originates from the thermal dissociation of 
stable molecules. 

We note that the key to understand the reason why the temperature dependence of Cy near 
slightly in the EEC regime is qualitatively different from the that in the unitary regime is the stabil¬ 
ity of preformed pairs. In the former EEC case where long-lived stable molecules (Nb - N/2) with 


10 
















a finite binding energy Ebind dominate over the system, thermal dissociation of molecules leads to 
exponential-like temperature dependence of various thermodynamic quantities. Because of this, 
the entropy S becomes a concave function of temperature, so that the specific heat Cy given by 
Eq. (fT3l) becomes an increasing function of the temperature. On the other hand, in the case of 
unitary regime with no stable molecule (A^b = 0)> because metastable preformed Cooper pairs are 
actually pairing fluctuations where formation and dissociation of preformed pairs repeatedly and 
frequently occur, the binding energy of such a fluctuating quasi-molecule is somehow ambiguous, 
especially when the molecular lifetime is very short. As a result, fluctuating metastable preformed 
pairs would not not give an exponential temperature dependence of S. However, as mentioned pre¬ 
viously, the growth of low-energy pairing fluctuations or metastable preformed pairs with nearly 
zero center of mass momentum near Tc decreases the entropy S. In addition, because this growth 
is more remarkable at lower temperatures near Tc, the entropy S becomes a convex function of 
temperature, so that Eq. (fT3]) gives the quite opposite temperature dependence of Cy to the case 
slightly in the BEC regime {(kyOsT^ 2*0.3) near Tc . 

The above-mentioned difference can be also understood on the viewpoint of the internal energy 
E. In the unitary regime, while the molecular picture is ambiguous, pairing fluctuations are known 
to induce particle-hole coupling Obll . leading to a pseudogap structure around the Eermi level. This 
phenomenon would lower the internal energy E (as in the case of the ordinary BCS state), which 
would become more remarkable at lower temperature near Tc, because of the enhancement of 
pairing fluctuations. As a result, E becomes a convex function of temperature, so that Eq. ([7]) gives 
the anomalous amplification of Cy in the unitary regime near Tc. On the other hand, slightly in the 
BEC regime where low-energy single-particle Eermi excitations near Tc are dominated by thermal 
dissociation of long-lived stable molecules, the internal energy would have an exponential-like 
temperature dependence, so that Eq. (|7]) gives the increase of Cy with increasing the temperature 


near Tc. 

Deep inside the BEC regime 2*0.7), we see in Eigs. |5];al)-(cl) that the enhancement 

of Cy is revived near Tc, although it is not so remarkable as the case of unitary regime. In the 
temperature region where Cy increases with decreasing the temperature, Eigs. [5];a2)-(c2) show that 
the system is dominated by long-lived stable pairs (Ab - A/2). In addition, Cy in this temperature 
region is well described by the specific heat C® of an ideal Bose gas with A/2 molecules, as shown 
in Eigs. [2al)-(cl). Thus, when one conveniently introduces the characteristic temperature T as 
the temperature at which Cv{T) takes a minimum value, the region Tc < T < f may be regarded 
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as an almost ideal Bose gas of long-lived stable molecules. 

We briefly note that such a Bose gas behavior of Cy can be also confirmed analytically. In the 
st rong coupling BEC regime, the particle-particle scattering matrix T(q, ivn) in Eq.® is reduced 


tom. 


Y{q,iVn) = 


Stt 


1 




(16) 


where /Jb = 2/u + In obtaining Eq. (fT^. we have used the well-known result in the BEC 
regime, // - = -IKlma^s^) «: -epllJ-lsl]. Substituting Eq. (fT^ into the internal energy 


in Eq.®, one has 


4m ' 4-m 


4m 


■ Fb 


F ^ 

^bind • 


(17) 


Here, we have assumed that all the Eermi atoms form tightly bound molecules, for simplicity. 
Since the specific heat Cy = {dFldT)yf^ is simply obtained from the first term in Eq. (fT71) . it is 
just the same as the specific heat in an ideal gas with N/I two-body bound molecules. 

With increasing temperature above f, the gradual decrease of the number Nb of stable pairs 
from A^/2, as shown in Eigs. [5ta2)-(c2), indicates the thermal dissociation of molecules. As shown 
in Eq. (fTSl) . this phenomenon naturally increases Cy, giving the deviation from C® seen in Eigs. 
[5l:al)-(cl). 

Plotting T in the phase digram of an ultracold Eermi gas in terms of the interaction strength 
and the temperature, we obtain Eig. As discussed above, this line physically gives the boundary 
between the region (NB) of an almost ideal Bose gas with A/2 non-condensed long-lived stable 
pairs and the so-called pseudogap regime (PG), where metastable preformed pairs dominate over 

3.8. We birefly point out that this 


the system. Particularly at T^, this boundary is at {kyasY^ - 

value is consistent with the previous result {k^asY^ - 0-75 [13911 . which was determined from the 
analyses of Eermi single-particle excitations. 

We note that the boundary between the pseudogap regime (PG) and the normal Bose gas regime 
(NB) has previously been given by T’ = 2|/i| (where /r < 0)[24, The background idea for 
this characteristic temperature is that 2|/i| eventually coincides with the binding energy = 
1 /(mol) of a two-body bound molecule in the BEC limit, so that stable molecules are expected 
to appear below T' ~ overwhelming thermal dissociation. However, comparing T' with T, 
one finds that they are actually very different, as T' f (although we do not explicitly show 
this comparison here). This indicates that, although stable pairs would start to appear around 


12 















FIG. 6: (Color online) Phase diagram of an ultracold Fermi gas with respect to the interaction strength 
{k^asY^ and the temperature T, scaled by the Fermi temperature Tp. The characteristic temperature f gives 
the boundary between the region (NB) of an almost ideal Bose gas with N12 non-condensed long-lived 
stable pairs and the pseudogap regime (PC), where the system is dominated by metastable preformed Cooper 
pairs or pairing fluctuations. T physically gives the boundary between the normal Fermi gas regime (NF), 
and PG. The region below is in the superfluid state. For comparison, we also plot the previous pseudogap 


temperature T* 112411 obtained from the density of states p(m), as well as the spin-gap temperature r J29tl 


determined from the spin susceptibility Xs- We note that Tc is only the phase transition temperature, f, T, 
T*, as well as Tg, are all characteristic temperatures, without being accompanied by any phase transition. 


T' ~ immediately mean the realization of a moleeular Bose gas. To obtain a gas 

of long-lived stable pairs, we need to further deerease the temperature down to f, at least on the 
viewpoint of the speeifie heat Cy. In this sense, the region between T' and t may be regarded 
as the erossover region between a gas of metastable quasi-moleeules and that of long-lived stable 
moleeules. 


We note that the physieal meaning of T is different from the previous pseudogap temperature 
r* [|24|] and the spin-gap temperature T, 1127112911 . because the latter two physically give the bound¬ 
ary between the normal Fermi gas regime (NF) and PG. In this regard, we point out that the specific 
heat Cy can also give the other characteristic temperature, which we denote f, corresponding to T* 
and Ts. As seen in Fig. |71 when one move to the weak-coupling BCS side ((kpasY^ ^ 0) from the 
unitary regime, the enhancement of Cy near T^. (which is caused by metastable preformed pairs) 
gradually disappears, and the temperature dependence of Cy is reduces to that in a free Fermi gas, 
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FIG. 7: (Color online) Calculated specific heat Cy in the BCS side ((A:Fai)“^ < 0) as a function of temper¬ 
ature T, scaled by the Fermi temperature Tp. Cy is the specific heat in a free Fermi gas in Eq. (fT^ . The 
characteristic temperature T is determined as the temperature at which Cy becomes minimal. 


given by 


4 = 2Xi- 


■P QJ ■ 


(18) 


Equation (fTSl) is proportional to T when T ^ Tp (see Eq. (fTTl) ). It approaches the classical 
Dulong-Petit law, Cy = 3Af/2 ll57n in the high temperature region. Thus, the temperature (= f) at 
which Cy takes a minimal value in Eig. |7]may be reasonably interpreted as the boundary between 
the normal Eermi gas regime (NE) and the pseudogap gap regime (PG) dominated by fluctuating 
metastable preformed Cooper pairs. Indeed, when we plot this characteristic temperature f in Pig. 
it is found to be consistent with T* and Eg, as expected. 


Pinally, we compare our result with the recent experiment on a ^Ei unitary Eermi gas 114 ill . Pig 


ure[8]shows that our result well explains the observed amplification of Cy near T^, indicating that 
the observed anomaly is due to metastable preformed Cooper pairs. However, Pig. [8] also shows 
that our result overestimates this enhancement near T^. In this regard, we recall that a finite spacial 
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FIG. 8: (Color online) Comparison of our result with the resent experiment on a ®Li unitary Fermi gas 14ill . 


In this figure, the temperature is normalized by T^. Cy is the specific heaf in a free Fermi gas. 


resolution inherent in this experiment in a trapped geometry eould lead to a possible suppression 
of the speeifie heat near r^ ll4ll] . We also point out that, sinee we deal with pairing fluetuations 
within the simplest NSR level, inelusion of higher-order strong-eoupling eorreetions beyond this 
approximation may also be important to eorreetly deseribe the behavior of the speeifie heat Cy, 
espeeially near Tc. Thus, we need further analyses to quantitatively explain this experimental 
result. 


IV. SUMMARY 

To summarise, we have diseussed the speeifie heat at eonstant volume Cy in the BCS-BEC 
erossover regime of an ultraeold eold Fermi gas. Ineluding pairing fluetuations within the frame¬ 
work of the strong-eoupling theory developed by Nozieres and Sehmitt-Rink, we elarified the 
temperature dependenee of this thermodynamie quantity over the entire BCS-BEC erossover re¬ 
gion above T^. In the unitary regime, we found that the speeifie heat is anomalously amplified near 
Tc, whieh is due to the appearanee of fluetuating metastable preformed Cooper pairs. Although 
this anomaly onee disappears as one goes to the BEC side, Cy was found to be again enhaneed 
near Tc with further inereasing the interaetion strength. We showed that this regime is dominated 
by long-lived stable pairs, and the enhaneement of Cy in this region agrees well with the ease of 
an ideal moleeular Bose gas. 
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Using these results, we determined the eharaeteristie temperature T whieh physieally distin¬ 
guish between the pseudogap regime where the system is dominated by metastable preformed 
Cooper pairs (or pairing fluctuations) and the region of an almost ideal Bose gas with Njl non- 
condensed stable pairs. 

From the temperature dependence of the specific heat in the BCS side, we also determined 
the other characteristic temperature f, which physically distinguish between the normal Fermi 
gas regime and the pseudogap regime. Using T and f, as well as the superfluid phase transition 
temperature T^, we obtained the phase diagram of an ultracold Fermi gas in terms of the interaction 
strength and the temperature, consisting of (1) the normal Fermi gas regime, (2) the pseudogap 
regime dominated by metastable preformed Cooper pairs or pairing fluctuations, (3) the region of 
an almost ideal Bose gas with N12 non-condensed long-lived stable pairs, and (4) the superfluid 
phase below Tc. Although t and f are not accompanied by any phase transition, they are still 
useful in considering the strong-coupling properties of an ultracold Fermi gas in the BCS-BEC 
crossover region. 

We note that, although the background physics of f is similar to that of the previous pseudogap 
temperature T*, which is determined from the density of states p(m), as well as that of the spin- 
gap temperature T^, which is determined from the spin susceptibility it is difficult to obtain 
the characteristic temperature corresponding to f from p(m) This is because they almost 

vanish in the strong-coupling BEC regime, due to the formation of tightly bound spin-singlet pairs 
with a large binding energy. In contrast, the specific heat is not suppressed in the BEC regime, 
so that we can safely determine t to identify the region consisting of stable pairs below T. In 
addition, the specific heat is known to exhibit singularity at T^. These advantages indicate that the 
specific heat is a useful quantity in constructing the phase diagram of an ultracold Eermi gas in the 
BCS-BEC crossover region. 

We note that we have included strong-coupling corrections within the simplest NSR theory in 
this paper. In this regard, while the NSR theory can describe the BCS-BEC crossover behavior of 
Tc, this strong-coupling theq^ is known to overestimate the pseudogap phenomenon associated 


with pairing fluctuations 1124. 


27h . Since the NSR specific heat at the unitarity overestimates the ob¬ 


served enhancement of Cy near T^. in a ^Ei Eermi gas (see Eig. [8]), a more sophisticated treatment of 
pairing fluctuations beyond the NSR theory would be necessary, in order to quantitatively explain 
this experiment. In addition, since the NSR theory completely ignores an effecvtive interaction 


between molecular bosons 1161 


6211 . it is also a crucial issue to clarify to what extent this molecular 
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interaction affects the characteristic temperature T (which physically gives the boundary between 
the region of (long-lived) stabe molecules and the region of metastable preformed pairs). For this 
problem, the so-called self-consistent T -matrix approximation ^ II would be useful. 

We also note that we have only dealt with the normal state above in this paper. Thus, ex¬ 
tension of the present theory to the superfluid phase below is also an interesting challenge. In 
addition, we have ignored effects of a harmonic trap in this paper. Although these effects should in 
principle be unimportant regarding the fact that the recent experimental result shown in Figj^rep- 
resents that of a uniform Fermi gas, it has been pointed by the authors of this experiment s ill that 
the trap geometry may induce an error in the temperature measurement. To quantitatively com¬ 
pare our result with the experiment data, we need to theoretically include this point. Although the 
pseudogap phenomenon associated with pairing fluctuations has recently attracted much attention 
in cold Fermi gas physics, the pseudogap temperature between the normal Fermi gas regime and 
the pseudogap regime has so far been mainly discussed. Thus, our results would contribute to the 
further understanding of BCS-BEC crossover physics in an ultracold Fermi gas on the viewpoint 
of the preformed pairs. 
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Appendix A: Expression for the number Ab of stable molecules 


To extract the contribution of stable molecules from the NSR term Ansr in Eq. Q, it is conve¬ 
nient to write it in the spectral representation, as 


■^NSR 



JmnB(m)pB(m). 


(Al) 


Here, = Yiq^B{q, oj) may be viewed as the molecular density of states, and the factor two 

means that each molecule consists of two Fermi atoms. The molecular spectral weight A^{oS) in 
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has the form, 


AB(q,oj) = --Im 
n 


r(q, m+) 


d 


d(2fi) 


U(q, aj+) 


(A2) 


In Eq. (IA2h . we have used the simplified notations, r(q, m+) = r{q,iv„ co + id), and 
n(q, CO+) = n(qr, iv„ to + id), where d is an infinitesimally small positive number. When the 
analytie eontinued partiele-partiele seattering matrix r(qr, to+) has a real pole at m = tOq, it ean be 
approximated to 

1 


r(q,m+) = = 


m 


Anas 




[a)+ - a)q]-£-Iliq,a)q) 


(A3) 


The eontribution of the pole at to = tOq to the number equation is evaluated by substituting Eq. 
(IA3I) into Eq. (IA2I) . Sinee the real pole tOq physieally deseribes the dispersion of a stable moleeule, 

^ ^n(qr, tOq) 

q:pole 




(A4) 


has the meaning of the number of stable pairs, where the summation is taken over real poles of 
r(q,to+). 

The eontribution N^c of seattering states to the number N of Eermi atoms is then given by 


A. 


Ansr - 2Ai 


B- 


(A5) 
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